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ABSTRACT 

We derive an analytic solution for the minimization problem in the geometric 
Baade-Wesselink method. This solution allows deriving the distance and mean 
radius of a pulsating star by fitting its velocity curve and angular diameter mea- 
sured interferometrically. The method also provide analytic solutions for the 
confidence levels of the best fit parameters, and accurate error estimates for the 
Baade-Wesselink solution. Special care is taken in the analysis of the various 
error sources in the final solution, among which the uncertainties due to the 
projection factor, the limb darkening and the velocity curve. We also discuss 
the importance of the phase shift between the stellar lightcurve and the velocity 
curve as a potential error source in the geometric Baade-Wesselink method. We 
finally discuss the case of the Classical Cepheid C Gem, applying our method to 
the measurements derived with the Palomar Testbed Interferometer. We show 
how a careful treatment of the measurement errors can be potentially used to 
discriminate between different models of limb darkening using interferometric 
techniques. 

Subject headings: Cepheids — stars: fundamental parameters — stars: oscilla- 
tions — techniques: interferometric — stars: individual {C, Gem) 



1. Introduction 



Since its conception, the Baade-Wesselink (BW) method (Baade 1926; Wesselink 1946) 
has been adopted as a preferred way to measure the distance of pulsating stars. In its 
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classical formulation, the distance modulus of a star is derived from the stellar physical 
radius, obtained by integrating the stellar radial velocity curve, and its effective temperature. 
The method has been applied to calibrate the Period-Luminosity relation for different classes 
of pulsators, among which RR Lyrae (McDonald 1977; Jones et al. 1992; Cacciari et al. 1992; 
Fernley 1994; Bono et al. 1994), Cepheids (see review by Gautschy 1997), Delta Scuti and 
SX Phoenix variables (Meylan et al. 1986). 

The main limitation of this technique, however, consists in finding a reliable observable 
yielding an accurate estimate of the stellar T^s (Bohm-Vitense et al. 1989). As a result, 
many variants of the BW method have been developed to circumvent this problem, by using 
different combinations of colors and band-passes (see e.g. Gieren et al. 1993; Laney & Stobie 
1995; Ripepi et al. 1997; Balog et al. 1997). Despite these ongoing observational efforts, 
discrepancies in the BW distances derived with different colors are still present (Laney & 
Stobie 1995; Krockenberger et al. 1997). 

Recent progress in interferometric techniques have made possible direct determination 
of angular diameter variations in nearby pulsating stars. This achievement has resulted in 
the development of a geometric version of the Baade-Wesselink method, which is in principle 
free of the limitations of the color-based techniques. Apphed to the galactic Cepheids, this 
method should eventually lead to a reliable zero point of the Cepheid distance scale (Sasselov 
& Karovska 1994). 

Attempts to measure the angular radial displacements of pulsating Cepheids were made 
using the IOTA interferometer (Kervella et al. 2001), the Palomar Testbed Interferometer 
(PTI, Lane et al. 2000, 2002) and the Naval Prototype Optical Interferometer (NPOI). The 
results of these observations show that indeed the BW method can be used to yield distances 
of nearby Cepheids with a much better accuracy than other methods, including geometric 
paraUax (Feast & Catchpole 1997; ESA 1997). 

The geometric BW method is limited in its accuracy by uncertainties due to observa- 
tional errors and model dependencies. These model dependencies derive from the need to 
translate the observed visibilities into angular diameters, which require accurate center-to- 
limb brightness profiles specific for the pulsating star. These profiles can be derived from 
hydrodynamic models of the stellar atmosphere (Marengo et al. 2002, hereafter Paper I). The 
conversion of radial velocities into radial displacement requires the knowledge of a projection 
factor which is computed by modeling line formation in the pulsating atmosphere (Sabbey 
et al. 1995). To these model-dependent uncertainties one should add the intrinsic errors 
of the interferometric measurements, and possible irregularities in the periodic pulsational 
behavior of the star. 
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In this paper we revisit the geometric BW method, in order to discuss the relative 
importance of various uncertainties, assessing their contribution to the accuracy of the final 
distance determination. We first derive an original analytical solution for the BW method 
fitting procedure, which simplify the determination of the best fit distance and average 
radius from the observational data, and allows a more transparent analysis of the individual 
error sources. We then apply our revised method on ( Gem PTI data (Lane et al. 2002), 
and we conclude with a discussion of the importance of the errors due to the pulsational 
model uncertainties (projection factor and limb darkening) and the phase shift between the 
lightcurve of the pulsator and the dynamical phases of the radial velocities. 



The geometric BW method allows deriving the distance and the mean radius of a pul- 
sating star by fitting the variations of its angular diameter measured interferometrically, 
with the radial displacement Ait!(0) of the stellar photosphere. The radial displacement is 
in turn derived by integrating the pulsational velocity over time: 



where Vj.{ct)) is the radial velocity, which should be corrected for the systemic velocity 7 and 
the projection factor p((p) to yield the pulsational velocity. 

As explained in Paper I the systemic velocity is calculated by requiring the conservation 
of the radius over one period. It takes into account the radial motion of the star relative to the 
Solar System, and physical inconsistencies in the measurement of the radial velocities with 
spectral data. An appropriate p-factor is computed by modeling spectral line formation for a 
given pulsating star, including spectral line asymmetries and the dynamics of the pulsation. 
Models for several classical Cepheids have been computed by Sabbey et al. (1995), from which 
the pulsational phase-dependent p-factors are derived. The variations of the radial velocity 
with the pulsational phase arc instead measured spectroscopically. A recent compilation of 
radial velocity data for 40 Cepheids have been published by Bcrsicr et al. (1994b), by using 
the CORAVEL spectrophotometer. The measurements arc provided as a Fourier expansion 
of Vr, in terms of the pulsational phase 0. The integration of the Vr{(f)) curve provides 
the required radial displacement Ai?(0), used in combination with the measured angular 
diameters in the geometric BW method. 

The BW distance D and average radius Rq of the pulsating star are solved as a two 
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Geometric Baade-Wesselink method 




(1) 



-4- 



parameter minimization: 



X^(i?o,^) = E 

i 

where the index i runs on the data points, Qi is the angular diameter measured at a certain 
phase 0j and Ai?j = Ai?(0j) is the radial displacement derived with equation 1 for the same 
phase. The terms o"j are the errors associated with each data point, including the uncertainty 
in the radial displacement. Note that when converting the time of each individual data point 
Qi into the pulsational phase 0j of the radial displacement curve, special care is required to 
take into account the phase shift observed between the stellar lightcurve and the radial 
velocity curve from which is derived. Note also that the diameters ^j, when measured 

with an interferometer, should consider the Limb Darkening (LD) of the stellar atmosphere. 
As in the case of the p-factor, this quantity is also affected by the stellar pulsation, and 
depends on the pulsational phase. In paper 1 a method was presented to compute wavelength 
and phase dependent LD, using the same hydrodynamic models from which the Sabbey et 
al. (1995) p-factor was derived. 

The fit for the geometric BW method is solved by minimizing equation 2 with respect 
to the average radius and the distance D. The minimum is found at the stationary 
point of X^(-Ro, -D): 



E 



(2) 



5roX' = 
dnx^ = 



2 _ . (3) 



which gives a system of two linear equations in D, quadratic in Rq: 



U - -^Ro + ^ 



BRo+E 



where the following coefficients are defined as: 
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The system in equation 4 describes two curves intersecting at the minimum . The 
angular coefficient of the fcst curve, which is a ffist order polynomial, is 2AjB. The second 
curve, when linearized to the ffi^st order, has an angular coefficient also very close to 2A/B. 
This means that the two curves described by the system in equation 4 are almost parallel. 
Their intersection will thus result in a large uncertainty along the D/Rq — 2A/B direction, 
and a much smaller error on the orthogonal direction. 

The system in equation 4 can be solved by direct substitution, obtaining first an explicit 
solution for Rq in terms of the coefficients defined in equation 5. When the solution for Rq 
is substituted back into the second equation, the quadratic terms cancel, leaving a unique 
solution for D. Thus the best fit distance D and mean radius are: 



The value of the for the best fit parameters is obtained by substituting Rq and D 
into equation 2, obtaining the following expression: 



The error region for the best fit solution is determined by the confidence levels of the 
X^(i?0) D). Assuming a gaussian error distribution, the 68% confidence level (giving the la- 
error) is defined by the equation: 




(6) 




Ro=Ro 



(7) 



J D=D 
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where Xggo/ = 2.33 for a two parameter fit (see e.g. Prodesen et al. 1979). The explicit form 
of the error contour in terms of Rq and D is obtained from equation 7 by substitution of 
Xmin with Xmin + ^ggo/ • Notc that the equation describes a tilted ellipse having the major 
axis oriented along D — Rq. 

Error bars for the Rq and D best fit solution, at a certain confidence level, can be 
obtained analytically by determining the intersection of the error ellipse with appropriate 
lines parallel to the Rq and D axis. Each line intersects the ellipse in two points (given by 
a quadratic equation). The intersections which provide the error bars are the ones where 
these two points coincide; there arc two such intersections for each axis, each of them given 
by the solution of a quadratic equation in D and Rq respectively. These solutions can be 
solved explicitly in terms of the coefficients defined in equation 5 as: 



^±^^^±±^W^E^ (9) 

ax 

where X is either D or Rq, and where the numeric coefficients in the two cases are: 




and 



2 {AE - CB) 
4 (Cf - AF) 



(10) 




S^_^G + ^(x- + x^8% 
BE + c{x- + xl^^yCG 
B' + f-(x^ + X^8%)-re 



(11) 



The best fit angular diameter is simply obtained as Oq = Rq/D. Note that the expres- 
sion for = Rq/ D is the equation for the major axis of the error ellipse. Lines parallel 
to the ellipse major axis are the regions of constant angular size for the generic Rq and D 
solutions; the error bars for ©o are thus given by the intersections of these lines with the 
error eUipse. The solution can then be found analytically, as in the case of the error bars of 
Rq and D, with the following coefficients: 
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3. Application to the Classical Cepheid C, Gem 

To illustrate the procedure described in the previous section, and analyze the relative 
importance of the various error sources in the final solution, we can consider the case of the 
Classical Cepheid ( Gem, recently observed with the PTl interferometer (Lane et al. 2002). 

The Uniform Disk (UD) angular diameters of^^'' obtained for ( Gem in the near-lR 
H-band using the PTl have been published in Table 4 of Lane et al. (2002). Each diameter 
results from a number of individual measurements (scans) performed at a certain date (in 
Julian Days, JD). The JD of each measurement is converted into the ( Gem lightcurve 

(i) 

phase (f))^ by using the published period P = 10.150079 days and the zero-phase epoch 
To = 2,444,932.736 (Lane et al. 2002). As mentioned previously, the fitting procedure 
requires reconciling the dynamical" phase (pv of the radial velocity curve with the lightcurve 
phase (pL of the observations. This is done by correcting for the phase shift A0 = (pL — (pv 
observed between the maximum luminosity and the minimum radius. Prom the data collected 
by Bersier et al. (1994a,b) we derive A0 ~ 0.28. This is the value we will use in the 
following discussion. Limitations of this approach for deriving the phase shift are discussed 
in section ssec-shift. 

To derive an accurate distance with the geometric BW method, the angular diameter 
measurements 9i should be corrected for the LD. Por clarity, however, we first solve the BW 
fit using the published UD diameter, deferring the discussion of the effect of the LD, and the 
errors associated with it, to a separate section. Por the same reason, we adopt the value of 
the p-factor of 1.43 used in Lane et al. (2002), which is the average value derived by Sabbey 
et al. (1995) over the ( Gem pulsational cycle, consistent with the model we adopted for 
computing the LD. 

As discussed before, the terms (Xj are the errors associated with each data point. They 
are the geometric sum of the individual errors for each observation: 



+ 



(AR/D) 



(13) 



where af ^ is the error of the interferometric measurement and al^^'^' is the error related 
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to the diameter variation. According to Lane et al. (2002) the errors in the interferometric 
PTI UD data vary between 0.01 and 0.06 mas (10-60% of the predicted amphtude of ( Gem 
pulsations, of ~ 0.1 mas), depending on the quahty of each measurement. Note that these 
errors refer only to the determination of the UD diameters. When computing the true 
distance of the star using the diameters corrected for the LD, we will also take into account 
the errors involved into the LD correction. 

An estimate of the error related to AR/D is more difficult to obtain. The ffi^st compo- 
nent in this error is due to the uncertainty in the radial velocity measurement. From Bcrsicr 
et al. (1994a) we get that the average error in the ( Gem velocity measurement is about 
0.4 km/s over velocities up to ~ 20 km/s, which takes into account possible irregularities in 
the pulsational behavior of the star. This is a relative error of up to ~ 2% (when the veloc- 
ity is maximum) , which translates to a maximum uncertainty in the expected stellar radius 
variation of ~ 0.002 mas for the predicted amplitude of the ~ 0.1 mas ( Gem pulsation. 
Another contribution to this error is the uncertainty in the p-factor; we use the same value 
of 1.43 ± 0.06 as Lane et al. (2002), which takes into account the variation of Sabbcy et al. 
(1995) p(0) with the pulsational phase. This adds a further ~ 4% relative uncertainty on 
the measured pulsation amplitude (equivalent to ~ 0.004 mas). The total error associated to 
the ( Gem angular diameter variation due to the uncertainties in the AR/D fitting function 
(again assuming that the two error sources are not correlated, and add as their geometric 
sum) is thus of ~ 4.5%. This is equivalent to a maximum error of ~ 0.005 mas for the 
predicted amplitude of the C Gem pulsations. Note that this error, in general, is smaller 
than the error of the current interferometric data, but may become important in the future, 
given the promising advances in the interferometric techniques. To include the af^^^^'' in 
the BW fit, we have first solved equation 6 in order to have a preliminary estimate of D, 
and then evaluated the a^-^^^^^ ~ (0.06 ■ ARi/D) for each point. Finally, we have repeated 
the fitting procedure considering both error sources. 

The best fit results are summarized in Figure 1, which shows the BW distance and 
average radius for the PTI UD data, and the error regions for the 90% (external oval) and 
68% (inner oval) confidence level. As expected, the error regions are ellipses with the major 
axis oriented along the Rq/D = const direction. The ellipses are very narrow, resulting in 
large errors, individually, in the best fit Rq and D, but a small error in their ratio Gq = Rq/D. 
This is not surprising, since the accuracy of the interferometric measurements is good enough 
to produce a reasonable fit to the average angular radius of the pulsating star; a much higher 
accuracy is instead required to measure the radial displacement, from which Rq and D are 
measured with the BW method. Equation 2 can in fact be written as: 
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Given that the maximum angular diameter variations are less than 10% of the average 
angular diameter, then for most phases Go ^ 2ARi/D. This means that the fit is essentially 
a one-parameter fit for ©o, which can be determined with good accuracy. The measurement 
errors plays a much larger effect in limiting the accuracy of Rq and D along Rq/ D — Qq. 

3.1. Phase Shift Determination 

Two basic questions remain open when discussing the accuracy of the fit using the 
BW method. The first concerns the assumption that the phases of the interferometric data 
and radial velocities are well synchronized, e.g. that there is not an unknown phase shift 
component between the luminosity lightcurve and radial pulsations. The second question 
concerns the LD correction, which we need to apply to the UD measurement to derive the 
correct BW distance. 

As mentioned before, the BW method requires defining a common phase reference for 
diameter observations and radial velocity data. This is usually provided by the lightcurve 
phases 0^, with the zero-phase set to coincide with the maximum luminosity. A complication 
can however arise when observable quantities arc related to phase dependent parameters 
which are derived from dynamic modeling of pulsating stars. Typical parameters are the 
p-factor (Sabbey et al. 1995) and the LD correction (Paper I). The traditional choice of 
the lightcurve phase is not convenient when dealing with time-dependent hydrodynamic 
models, e.g. the ones from which the ]9-factor and the LD are computed. The emergent 
flux in a given waveband is in fact a derived quantity (which can only be computed with 
extensive radiative transfer modeling) , and is usually not available in hydrodynamic numeric 
computations. A better approach in this case is to choose the zero-phase reference based on 
hydrodynamic quantities directly tied with observables, as the radial velocity Vr- This leads 
to the "dynamical phase" 0y mentioned before. 

The complication in reconciling the lightcurve and dynamical phase references is due to 
the phase shifts occurring between the zero phases in the two systems. This quantity can be 
obtained by carefully measuring the lightcurve and velocity curve of the pulsating star, as 
done by Bersier et al. (1994a,b). This procedure yields for ( Gem the quantity ~ 0.28 
used in the previous section; different values should be expected for different stars, since the 
phase shift is determined by the unique dynamics of the pulsator. A different approach was 
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foUowed by Lane et al. (2002), in which the phase shift was derived as a separate quantity 
in the BW fit, on the assumption that having a third free parameter would produce a better 
agreement between the measured angular diameters and the radial displacement curve. 

To assess the effects on the fit results by an independent determination of the phase 
shift, we repeated the BW fit of the ( Gem PTI UD diameters leaving this parameter free. 
The results are showed in Figure 2, where the best fit is plotted in the case of a free phase 
shift (dashed line), and a fixed phase shift (solid line) as determined from Bersier et al. 
(1994a,b) data. The best fit results with their errors are shown in Table 1. 

The results show that the small difference in the phase shift, although "visually" im- 
proving the fit of the data points, does not change significantly the best fit results. The 
difference in i?o, D and ©q is well below the error bars resulting from the fit, which suggests 
that the shift itself can be a statistical deviation due to the uncertainties of the data points. 
It is however interesting to investigate the possible causes of this discrepancy, since an unac- 
counted phase shift would have profound implications in terms of the evolution of pulsating 
stars. 

Since the actual phase of interferometric data acquired at certain JD is determined 
by folding many pulsational periods starting from the zero phase epoch, an uncertainty in 
the period determination can play a significant role in changing the measured phase shift. 
The typical uncertainty for the Cepheid periods tabulated in Bersier et al. (1994a) is of the 
order of ~ 10~^ days. The zero phase epoch for ( Gem is around day 2,444,933, while the 
observations were made between JD 2,451,605 and JD 2,451,896. Clearly the 300 days span 
in the time during which the observations were collected is not important since it corresponds 
to less than 30 periods. The time lag between the zero epoch and the data acquisition is 
instead of ~ 7000 days, which is approximately 690 periods. Given the accuracy of the 
( Gem period, the possible shift is of the order of ~ 0.07 days. This is less than a 0.01 phase 
shift error, which is one order of magnitude less than the unaccounted phase shift derived 
by fitting the PTI data. Therefore, if the measured phase shift is real, it is not explained by 
the current uncertainty in the pulsational period. 

Berdnikov et al. (2000) analyzed archival data of the maximum luminosity epochs for 
a number of Classical Cepheids, in search for secular variations of their period due to evo- 
lutionary changes. In the case of C Gem, they found a decrease of 0.07 days over a time of 
about 5000 pulsational cycles. This means a decrease of the period of 1.4 • 10~^ days per 
period, which is much less than the period uncertainty adopted in this paper and in Lane et 
al. (2002). 

Another unaccounted source of possible variations in the phase shift involves the possi- 
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bility that the interferometric angular diameter may be shghtly different (and off-phase) from 
the radius obtained by integrating the radial velocity measured from spectral lines. This dis- 
crepancy is mostly taken into account by the p-factor, but there may be unaccounted effects 
introducing an extra phase lag. Finally, even though the period may be stable for better 
than 10^^ days per cycle, there may still be a variable phase shift between maximum light 
and minimum radius, causing a drift of the AR{(pv)/D curve with respect to the lightcurve 
phase (pL. This is a possibility that can potentially affect the phase dependent quantities 
involved into the calculations using the BW method, among which the projection factor 
and the LD correction derived from hydrodynamic models, which require the exact knowl- 
edge of the phase relation between models (computed in terms of the dynamical phase) and 
observations (dependent on the lightcurve phase). There are no indications, however, that 
such effect is present. Our best fit results in Table 1, in fact, suggest that the uncertainty in 
the phase shift only affects the best fit diameter for less than 0.001 mas, which is less than 
1% of the estimated amplitude of the ( Gem pulsation. This error is small enough to be 
ignored when applying the BW method, compared to the other larger uncertainties in the 
data. 

These considerations suggest that, given the current accuracy of interferometric data, 
we are not in the position to provide an independent measurement of the phase shift with 
better quahty than the available data from Bersier et al. (1994a,b). The effect measured 
with the independent fit is very small and does not affect significantly the fit of Rq and 
D. For these reasons we conclude that a two-parameter fit for the BW method is currently 
preferable than leaving the phase shift as a third free parameter. 

However, we note that possible variations in the maximum luminosity - minimum ra- 
dius phase shift should be taken into consideration when the precision of interferometric 
measurements, and the accumulation of good quality data over long period of time, make 
their direct observation feasible. Even though unpredicted phase shifts may appear unlikely 
with our current data, their possible occurrence should be monitored, as they can provide 
important insights on the mechanics of pulsations, and their relation for the stellar evolution. 

3.2. Limb Darkening 

The remaining question concerning the accuracy of the BW fit involves the limb darken- 
ing. The LD correction is necessary in order to take into account the non uniform brightness 
of the stellar disk, resulting from the different depths in the stellar atmosphere probed by 
different lines of sight. When a star is partially resolved by interferometric observations, 
the existence of LD induces a change in the measured visibility with respect to the simpler 
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case of a uniform brightness disk. The translation of interferometric visibihties into angular 
diameters takes into account this effect, by fitting the data with limb darkened model visi- 
bilities. Alternatively, when the original data have already been fitted with a UD model, a 
LD correction can be applied. 

As in the case of the p-factor, the LD corrections can be obtained by solving the dy- 
namic structure of the pulsating stellar atmosphere, and then computing a complete radiative 
transfer model to derive the center-to-limb intensity profile. Using this approach we have 
derived a method to compute LD profiles for pulsating Cepheids, presented in Paper 1. We 
have then computed specific phase- and wavelength-dependent LD corrections for ( Gem, 
taking into account second-order accurate one- dimensional hydrodynamic calculations per- 
formed in spherical geometry, and a full set of atomic and molecular opacities (Marengo et 
al. 2003, hereafter Paper II). 

To convert the UD diameters Of^^"* used in the previous section into the LD ones, one 
has to divide by the LD correction fc(A, 0): 



niUD) 

Figure 3 shows the LD corrections we presented in Paper 11 for ( Gem. The top curve 
shows A; as a function of the pulsational phase for the PTl H-band. The model LD correction 
does indeed change with phase, with a ±0.3% variation around the average value k = 0.979. 
The largest change occurs close to minimum radius, when a shock wave crosses the Cepheid 
atmosphere (see Paper I). Barring systematic errors in our models, we can assume the 
uncertainty in the phase dependence as the total error on our LD correction estimate. Prom 
Paper II we have that this uncertainty results in an extra error of cr^^^^ — 0.003 mas (~ 3% 
of the ( Gem pulsational amplitude), which is small with respect to the interferometric data 
errors. The BW distance computed with this correction is shown in Figure 4 and in Table 2. 
For comparison we have also solved the ( Gem BW fit using the LD correction k 0.96±0.01 
(dashed line in Figure 3), derived from tabulated values computed by Claret et al. (1995) for 
static yellow supergiants. This values have been used by Lane et al. (2002) to derive the PTI 
BW distance for ( Gem published in their paper. Given the quoted uncertainty, this results 
in an extra error source of cr-^^^ ~ 0.02 mas (as much as 20% of the pulsational amplitude). 

Figure 4 shows that, within the quoted error uncertainties, the two solutions for the BW 
fit with different LD correction are mutually exclusive, as the respective la error regions 
do not intersect. In fact, the two LD solutions are also separate from the BW distance 
computed for UD diameters. This shows that, even if the final error range in the best fit 
i?o and D includes all three solutions, interferometric techniques can potentially be used to 
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discriminate between different models of limb darkening. Alternatively, this result can be 
seen as an indication of how an independent determination of the LD will result in a much 
better measurement of pulsating star distance with the BW method. 

4. Discussion and conclusions 

Table 3 summarizes the error balance of the geometric BW method. For each error 
source, the table provides the resulting uncertainty in the determination of the Cepheid 
angular diameter, in mas and as a fraction of the ( Gem 0.1 mas pulsation amplitude. The 
last column of the table shows how these uncertainties affect the final determination of the 
BW distance. The statistical error in the interferometric data is still the main source of 
error in the method, leading to the final la error bars of ~ 14% in our fit. The systematic 
errors (related to the radial velocity and phase shift measurements, the p-factor and the limb 
darkening) still play a secondary role. Note that the effect of the relative uncertainty of the 
systematic error sources, which we have estimated for the amplitude of the stellar pulsation, 
should be scaled down by roughly one order of magnitude to estimate their effect on the 
best fit distance, since the BW distance depends on the stellar radius which is more than 10 
times larger than the pulsation amplitude. 

The second error source in importance is the determination of the average value of the 
LD correction. The difference between the best fit value of the BW distance obtained with 
our ( Gem model and with the LD coefficient derived by Claret ct al. (1995) is ~ 3%. Our 
( Gem fit shows that using a LD correction that is not appropriate for the observed Cepheid 
can introduce an error as large as when not using LD at all (e.g. using a UD model). This 
stresses the importance of reliable modeling of the Cepheid atmosphere even when using 
interferometric data with the currently best available quality. 

The table also shows that the uncertainties in the radial velocity, phase shift, p-factor 
and phase dependence of the LD correction at near-IR wavelength, all together are respon- 
sible for less than a 1% error in the BW distance determination. These uncertainties plays 
only a secondary role compared to the measurement errors. However, with the constant 
progress in the development of the interferometric techniques, the visibility errors are bound 
to be reduced, and at that point these other factors, as discussed above, will become essential 
to fulfill the potential of the BW method. This will be especially true for interferometers 
operating at visible wavelengths. Figure 3 shows the LD correction we have computed for 
C, Gem with our model in the V-band (bottom solid line). The plot suggests that at certain 
phases in which the hydrodynamic effects play an important role (close to minimum radius 
and in presence of Shockwaves) the interferometric measurement of the Cepheid pulsation 
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amplitude can be affected by as much as 20%. This may result in a further uncertainty of 
up to 2% in the BW distance determination, which can only be corrected with an accurate 
time-dependent hydrodynamic modeling of the Cepheid atmosphere. 

Even considering all the error sources discussed in the previous sections, the current 
error bar for the C Gem BW diameter is already two times better than the quoted error 
in parallax measurements for this star, which is ~30% (Hipparcos, ESA 1997). A careful 
analysis of the error regions in the best fit plane, shows that the error in the best fit angular 
diameter is in fact much smaller than the individual errors in Rq and D, leaving room for 
dramatic improvements in the distance measurement if independent constraints arc set on the 
stellar radius. These can be derived by modeling the structure of the pulsating atmosphere. 

The results presented in this work show that a better determination of the limb darken- 
ing of the pulsating star can already lead to different estimates of the distance and average 
radius which are mutually exclusive. This justifies the recent efforts in producing accurate 
predictions for the model dependent quantities needed by the BW method, by means of de- 
tailed self-consistent hydrodynamic models of the pulsating atmosphere (see e.g. Paper I and 
II). An independent determination of the LD correction for nearby pulsating stars, which 
will be attainable once interferometers with baselines of several hundred meters will become 
operative, will thus provide a direct test for such models. This will open the road for a large 
scale application of the geometric BW method, to derive the distances of a large sample of 
pulsating stars, and thus attain the goal of an accurate calibration of their PL relation. 
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TABLE 1 

BW BEST FIT OF C GEM PTI H-BAND UD DATA 



fixed A(/) free 



A0 


0.280 


0.294 


Ro [R©] 


66 1+8-6 


66.21^1 


D [pc] 


3751^^ 


376^^^ 


Go [mas] 


1.630±0.007 


1.629±0.007 
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TABLE 2 

BW BEST FIT OF C GEM PTI H-BAND LD DATA 





k = 0.979 


k = 0.96 


Rq [Rq] 

D[pc] 
©0 [mas] 


67.0+19^ 
3721^^ 
1.665±0.007 


66.5^^-^7 
362i^^ 
1.699±0.007 
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TABLE 3 

INDIVIDUAL ERROR CONTRIBUTIONS SUMMARY 



Error source 


<T [mas] 


% of C Gem ampl. 


% of BW dist. 


6i measurement 


0.01 - 0.06 


10 - 60% 


14% 


Radial velocity 


0.002 


2% 


0.2% 


p-fact 


0.004 


4% 


0.4% 


Phase shift 


< 0.001 


< 1% 


< 0.1% 


LD average value 


0.02 


20% 


2% 


LD phase variations (H band) 


< 0.003 


<3% 


< 0.3% 


LD phase variations (B band) 


< 0.02 


< 20% 


<2% 
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Fig. 1. — BW fit of C Gem distance and mean radius from PTI H-band UD data. Inner 
contour is 68% confidence level (Icr), and outer contour is 90% confidence level error. 

Fig. 2.— Best fit UD diameter of C Gem PTI H-band data from Lane et al. (2002). The 
solid fine is the fixed phase shift as determined by Bersier et al. (1994a,b). The dashed fine 
is instead the best fit with the free phase shift. 

Fig. 3. — Phase dependent LD corrections for ( Gem computed by the model presented in 
Paper II. The LD correction is shown for the PTI H-band in the near-IR, and for the optical 
V-band. 

Fig. 4. — Best fit parameters and error regions for ( Gem PTI H-band data. Top curve is 
the UD data, middle curve is our best fit for model LD data and the bottom curve is the 
result obtained using a fixed LD correction of /c ~ 0.96 as in Lane et al. (2002). The inner 
region is the 68% confidence level of the fit (Icr), while the outer is the 90% confidence level. 
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